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Q I We present an investigation of the ultraviolet and X-ray spectra of the Seyfert 1.5 

galaxy Markarian 817. The ultraviolet analysis includes two recent observations taken 
with the Cosmic Origins Spectrograph in August and December 2009, as well as archival 
■ spectra from the International Ultraviolet Explorer and the Hubble Space Telescope. 

Twelve Lya absorption features are detected in the 1997 GHRS and 2009 COS spectra 
- of these, four are associated with high-velocity clouds in the interstellar medium, 
four are at low-significance, and the remaining four are intrinsic features, which vary 
between the GHRS and COS observations. The strongest intrinsic absorber in the 1997 
spectrum has a systemic velocity of ~ — 4250kms~^. The corresponding feature in the 



o 

5-H ■ COS data is five times weaker than the GHRS absorber. The three additional weak 

^ ■ (equivalent width from 13-54 mA) intrinsic Lya absorbers are at systemic velocities of 

— 4100kms~^, — 3550kms~^, and — 2600kms~^. However, intrinsic absorption troughs 
from highly ionized C IV and N V, are not detected in the COS observations. No ionized 
I absorption signatures are detected in the ~ 14 ks XMM-Newton EPIC spectra. The 

Q>^ I factor of five change in the intrinsic Lya absorber is most likely due to bulk motions 

I in the absorber, since there is no drastic change in the UV luminosity of the source 

' from the GHRS to the COS observations. In a study of variability of Mrk 817, we 

find that the X-ray luminosity varies by a factor of ^ 40 over 20 years, while the UV 
. continuum/emission lines vary by at most a factor of 2.3 over 30 years. The variability 

of the X-ray luminosity is strongly correlated with the X-ray power-law index, but no 
correlation is found with the simultaneous optical/UV photometry. 
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Introduction 



Long-term and often short-term variability in both flux and spectral shape, across the electro- 
magnetic spectrum, is a fundamental characteristic of the spectra of active galactic nuclei (AGN). 
While variability is ubiquitous in AGN, our understanding of long-term spectral changes in local 
AGNs is the result of studies of a few dozens of sources focusing primarily o n variations in broad 
optical emission lines in conjunction with reverberation mapping studies (e.g. JPeterson et al.ll2004l ) 
for local AGNs {z < 0.1). It is important to extend these studies in temporal space (e.g., acquiring 
a larger sample of observations over longer periods of time), the wavebands covered (e.g., multi- 
wavelength studies), and the number /type of sources sampled, in order to understand the physics 
behind the AGN emission - the accretion processes, nature of outflows, and associated feedback. 

In an effort to extend the variability studies of individual AGN, we present an analysis of the 
available UV and X-ray observ ations of Mrk 817, a nearby Seyfert 1.5 galaxy with a redshift of 
z = 0.031158 (|Falco et al.lll999l ). Its optical spectrum exhibits complex emission in the narrow-line 
region (NLR) and broad-line region (BLR), with outflowing NLR [O III] AA4959,5007A emission 



measured at a systemic velocity of —450 km s ^ ( Die et al.ll2006l ). In the UV, observations with 



both the Hubble Space Telescope (HST) and the Far Ultraviolet Spectroscopic Explorer (FUSE) 
also revealed outflowing gas. HST Goddard High- Resolution Spectrog raph (GHRS) obse rvations 
from 1997 of Lya indicated the presence of an intrinsic AGN outflow (jPenton et al.ll200d ). Later 
observations with FUSE, in 2000-2001, show weak absorption in O VI with a radial systemi c 
velocity component of — 4100kms~^, the largest offset measured in a Seyfert ( Dunn et al.l 120071 ). 
Analysis of multiple observations from FUSE shows that the weak absorption is variable in tha t 
the equivalent width of the absorber increases with decreasing continuum flux (jPunn et al.ll2008l ). 
Additionally, a second and third absorption component {v ~ —3000 kms~^) appeared in the final 
FUSE observation, potentially due to a change in ionization. 

Owing to its UV brightness and s i ght line through a complex of Galactic high velocity clouds 
(HVC Complex C. iGibson et ahlboOll : IColhns et ahlbood ). Mrk 817 was a selected target for the 
guaranteed time observations with the HST Cosmic Origins Spectrograph (COS). Two high signal- 
to-noise ratio observations were taken with COS, separated by ~ 5 months. Additionally, UV 
observations of this source taken with the International Ultraviolet Explorer (lUE) and HST GHRS 
allow a search for variability over a period of nearly 30 years. Details of the UV observations, are 
presented in Section [2j We include an analysis of the intrinsic AGN properties, including the UV 
continuum and emission-line properties from each of these observations. In Section [3l we present 
an analysis of X-ray spectra of Mrk 817, along with simultaneous optical/UV photometry, from 
XMM-Newton and Swift. Although the source is bright in the X-ray band, it has only recently been 
targeted for observation in the X-ray with XMM-Newton. Finally, we present our conclusions on 
the spectral analysis and UV and X-ray variability of Mrk 817 in Section [H We assume a standard 
cosmology with Hq = 75kms~^ Mpc~^, throughout. 
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2. The Ultraviolet Spectra 

In this section, we include details of our analysis of the available UV spectra of Mrk 817. In 
S 12.11 we describe the observations from COS, GHRS, and lUE. An analysis of the continuum and 
emission line properties from the COS observations is included in § 12.21 In § 12. 3| we compare the 
COS spectra with the GHRS and lUE spectra - including detected Lya absorbers in § 12.3.11 and 
the continuum/emission-line properties in § 12.3.21 



2.1. Observations 



Mrk 817 was observed twice with HST/COS with both the G130M (1135 A < A < 1480 A) and 
G160M (1400 A < A < 1795 A) medium-resolution {R = A/AA 18,000) gratings. The first epoch 
of observations (2009, August 4) were part of the Early Release Observations (ERO) program (PI: 
Noll), with four exposures in each grating for a total exposure time of 2162 s (G130M) and 1410s 
(G160M). The second set of exposures (2009, December 28) were taken as part of the COS Guaran- 
teed Time Observations (GTO) program (PI: Green) and totaled 1264 s and 1600 s in G130M and 
G160M, respectively. Both sets of observations used four different central wavelength settings per 
grating to dither instr umental features a n d provide continuous sp ectral coverage across the entire 
COS FUV band. See ICreen et al.l (|201ol ). lOsterman et"al] (|201ol l. and the COS HST Instrument 



Handbook (jPixon et al.ll20ld ) for details on the COS instrument and on-orbit performance. 



The ERO data were obtained during the Servicing Mission Orbital Verification (SMOV) period 
before the instrument reached final focus. However, the pre-focus data are indistinguishable from 
that taken after correct focus had been achieved, except in the case of the very narrowest of 
absorption features (AA < 0.1 A), none of which are analyzed in this work. Another symptom 
unique to the early SMOV data was the "divot and clod" features which occurred near the blue end 
of each detector segment; charge resulting from photons hitting one part of the detector appeared 
in a different part. The high voltage on the detectors was reduced on 2009, August 17, and the 
instrumental feature disappeared. For the Mrk 817 ERO observations, the divot and clod were 
carefully fitted with Gaussian profiles and normalized to the fiux level in adjacent portions of the 
spectrum. The wavelength region affected by this feature and its subsequent removal is relatively 
small, and much smaller in both amplitude and width than the AGN emission features analyzed 
here. Thus, their effect on our scientific conclusions is minimal. 

All data were reduced with CALCOS v2.11f. Partial fiat-fielding, alignment, and coaddition 
of the processed exposures were carried out using IDL routines developed by the COS GTO team 
specifically for COS FUV datfl Briefly, each exposure was corrected for narrow, ~ 15%-opaque 



^IDL routines available at jhttp: //casa. Colorado . edu/$\sim$da]if orth/costools .html| and described in 
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shadows from ion repeller grid wires. The local exposure time in these regions was reduced to give 
them less weight in the final co-addition. Similarly, exposure times for data at the edges of the 
detectors was decreased to de- weight these data in the final co-addition. With four different central 
wavelength settings per grating, any residual instrumental artifacts from grid-wire shadows and 
detector boundaries should have negligible effect on the final spectrum. 

Next, strong interstellar features in each exposure were aligned via cross-correlation and in- 
terpolated onto a common wavelength scale. The wavelength shifts were typically on the order of 
a resolution element (~ 0.07 A, ~ 17 kms^^) or less. The coadded fiux at each wavelength was 
taken to be the exposure-weighted mean of fiux in each exposure. 

The exposures from the ERO and GTO epochs were coadded separately and are presented in 
Figure [H To quantify the quality of the combined data, we identify line-free continuum regions, 
smooth the data by the seven-pixel resolution element, and define S/N = mean(flux)/stddev(flux). 
In the 1150-1200 A band, we measure S/N ~ 37 per resolution element in the ERO and S/N ~ 26 
per resolution element in the GTO observation. Similarly, in the 1500— 1550 A band, we measure 
S/N ~ 32 per resolution element in the ERO and S/N ~ 29 per resolution element in the GTO 
observation. 

In addition to the COS spectra, we also include archived lUE data from the Multimission 
Archive at STScI. We downloaded large-aperture, low-dispersion grating observations of Mrk817 
for two observations separated by approximately 8 months. The observation IDs are SWP 15447 
and SWP 17442 with observation dates of 1981 November 7 and 1982 July 18, respectively. The 
exposure time for each of these observations is 16200 s. The archived data were downloaded in 
the finished processed form. These spectra cover the AGN rest-frame wavelength range from 
~ 1120 - 2010 A. 

Finally, we include the HST GHRS spectrum of Mrk817. The UV spectrum was observed on 
1997 January 12 under project ID 6593 (PI: Stocke). Three exposures were taken with the G160M 
grating for a total exposure time of 25 ks. The wavelength coverage of the GHRS spectrum is only 
36A, cover i ng th e intrinsic Lya emission region. We use the processed spectrum presented in 



Penton et al.l (|200d ). which includes the details of the calibration and spectral processing. 



We summarize the details of all of the UV observations in Table [U including the observation, 
observation date, exposure time, spectral coverage, and spectral resolution for the COS, GHRS, and 
lUE spectra. Additionally, Figure [1] illustrates all of the UV spectra by order of their observation 
date. 



Danforth et"ai] (|201ol 'l 
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2.2. Analysis of the COS Spectra 

For the analysis of each o f the COS observations of M rk 817, we first corrected the spectra for 
Galactic reddening, using the iDickev &: LockmanI (1990) value of Nhi = 1-15 x lO^'^cm"^, which 
corresponds to E(B — V) = 0.022 following the relationship between hydrogen column density and 
interstellar reddening derived b y iPredehl &: SchmittI (Il995l ). We used this reddening value with 



the selective extinction curve of lCardelh et alJ Hwm and the near-UV correction fromlO'Donnel 



(j 19941 ) and then shifted the spectra into the AGN rest-frame. In Figure [H we show the resultant de- 
reddened COS spectra (along with the de-reddened GHRS and lUE spectra). The most prominent 
features are intrinsic AGN (and geo-coronal Lya) emission lines, and no strong intrinsic absorption 
is detected embedded in the broad emission line features. There are, however, a number of weak 
absorption lines in Lya, which we will discuss in 



We next chose continuum regions free of emission/absorption features to fit a power-law con- 
tinuum of the form: Fx{X) = ao x (A/1750 A)~", where oq is the flux density in units of erg s^^ 
cm "2 A"^ at 1750 A and a is the power-law index. Our choice of 1750A was arbitrary - correspond- 
ing to a wavelength free of prominent emission/absorption at the end of the COS bandpass. For 
our continuum regions, we chose 6 regions free of emission/absorption features from 1120-1700A, 
with individual region sizes from 10-50A. The same continuum regions were used for both the ERO 
and GTO spectrum. The resultant best-fit parameters, for both continuum and emission line fits, 
were determined using ^^-minimization along with the pyminuitj^ python interface to the MINUIT 
- Function Minimization and Error Analysis cod^ . Results of these fits are shown in Table [21 
with recorded errors corresponding to the l-cr level. We find that the power-law index is unchanged 
between the observations (a ~ 1.30 ±0.05). This slope is consistent with values obtained for other 
local AGN; for i nstance, the average slope fou nd from an analysis of HST obse rvations of narrow - 
line Seyfert Is bv IConstantin fc Shieldj ^200^ is a or Fa oc A^^^o and lScott etall (|20o3) 

find the composite EUV spectrum from FUSE of z < 0.1 sources as F,^ o c or F \ oc X~^'^^ 



We note that the slopes of higher redshift AGN tend to be different (i.e., iTelfer et al.l (120021 ) find 
-Fi/ oc u~^''^^ or Fx oc A"*^'^^ from the HST spectra of 184 quasars with z > 0.33). While the slope 
did not change appreciably between the COS ERO and GTO observations, the continuum level did 
vary, with the GTO observation being 30% lower than the ERO observation. 

The strongest features in the spectra of Mrk 817 are intrinsic AGN emission lines of Lya, the 
N V doublet, the Si IV doublet (blended with O IV]), and the C IV doublet. We fit the emission 
lines in these observations with a combination of Gaussians. We chose Gaussians as the simplest 
mathematical representation of the lines and use them for a basic characterization of the emission- 
line properties. Since the C IV doublet (AA1548.40, 1550. 77A) is the strongest feature that is not 



^http: //cod e .google . com/p/pyminuit7] 



^Copyright CERN, Geneva 1994-1998, with the 
[http: //wwwasdoc .web . cern. ch/wwwasdoc/minuit/minmain.html. 



MINUIT 



manual 



available 



at 
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blended with nearby emission lines, we first fit C IV with a broad doublet and intermediate doublet 
component (four Gaussians). Since the ratio of the statistical weights of the blue to red line is 
2 : 1, we fixed the flux of the blue line to be twice that of the flux in the red line. We also fixed 
the separation of the lines in the doublet to the laboratory measured value of 2.6 A. Best-fit values 
for the C IV doublet, including the velocity offset, velocity width, equivalent width, and line flux, 
as well as for the other strong emission lines described below, are shown in Tables [3] and [H 

To fit the Si IV doublet (AA1393.76, 1402.77 A), we fixed the velocity width of the hue to that 
found in the C IV fits. As for the C IV doublet fits, we fit the Si IV emission region with four 
Gaussians - a broad and intermediate component for each line in the doublet - fixing the difference 
between the line centers of each doublet component (broad or intermediate) to the laboratory 
measured difference between the lines in the doublet (9.015A) and the flux of the blue line to be 
twice that of the red line (the ratio of the line strengths is 2 : 1). We did not attempt to fit the 
Si IV and O IV] emission separately. 

Since Lya (A1215.67A) and the N V doublet (AA1238. 82, 1242.80) are blended, we fit these 
emission features simultaneously. We first masked strong Galactic ISM lines in this region, including 
the absorption features from N V, S II, and Si II. We fit the N V doublet with a broad and 
intermediate doublet composed of four Gaussians, fixing the width of the lines to be the same as 
that for the C IV doublet. As for C IV and Si IV, we fixed the separation of the lines in the doublet 
to the laboratory measured value (3.980 A) and the blue line of the doublet to twice that of the 
red line (the ratio of the line strengths is also 2 : 1 for N V). The asymmetric shape of the Lya 
emission required us to fit this feature with a combination of three Gaussians. We fixed the full 
width half- maximum for two of these Gaussians to the values obtained from the C IV fits. 

We plot the best-fit emission line models for the ERO and GTO spectra in Figures [2] and [3l In 
general, we find that our approximations to the emission lines with Gaussian models are quite good. 
The greatest discrepancies between the model and emission line-region occur, not surprisingly, in 
the Lya profile. For both the ERO and GTO spectra, we find that the peak of the line emission is 
not fully accounted for with our model. This is due to the complications in fitting the peak emission, 
which is coincident with strong Galactic S II absorption (AAA1250.58, 1253.81, 1259.52, 1260.422 A) 
that had been masked from the data for the spectral fitting. 

Comparing the results of our emission line fits between the ERO and GTO spectra, we find 
that the measured full- width-half- maximum (FWHM) values are similar for the C IV doublet. We 
find a best-fit for both spectra with velocities of the two Gaussian components corresponding to 
~ 3000 km s~^ and ~ 1000 km s~^. This is consistent with C IV FWHM values from the Palomar 



sample, where iBaskin &: Laon ( 20051 ) find a range of values from 1600-10250 kms ^ for their sample 



of 87 AGN at z < 0.5 AGN. 

The measured velocity offsets for each of the measured line profiles do not agree, however, 
either between observations or different lines in the same observation. For C IV, we find shifts in 
the line profiles from ~ —50 kms~^ to ~ —340 kms~^. These values are in line with what is found 
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bv IShang et alj ( 2007 ) (who find a mea n shift of —445 ±513 kms~^) in the Palomar- Green X-ray 
sample. As found in lShang et alJ (|2007l ). we find that C IV components with larger blue-shifts also 
have lower EW (i.e., in component 1, the GTO observation has a larger EW and a smaller velocity 
shift than component 1 of the ERO observation). 

We find that the Si IV emission is red-shifted by ~ 200 - 800 kms"\ while the N V emission 
is slightly red-shifted to blue-shifted by ~ 20 to —220 kms~^. We expect error in the Si IV 
measurements, since we did not de-blend the Si IV and O IV] emission. The N V measurements 
suggest a smaller shift than the measured C IV ofi^-sets. 

Velocity shifts in Lya indicate that the strongest component (1) is blue-shifted by ~ —1000 kms~^, 
with an additional red-shifted (800 kms~^) component (2). A third component, whose FWHM 
and EW changed slightly between observations {AFWHM ~ 100 kms^^ and AEW O.sA), 
has a blue-shift of — 556kms~^ in the ERO observation and — 57kms~^ in the GTO observa- 
tion. Typical off-sets betw een Lya and C IV (fLya — ^civ) for radio-selected quasars are from 
70-350 kms~^( Wilkeslll986l ). While we find much larger velocity shifts between Lya and C IV, it 
is difficult, however, to make a direct comparison since our results are based on fits using multiple 
Gaussians as opposed to single Gaussian fits for each emission line. 

Additionally, we find measurable diff^erences in the emission line flux between the ERO and 
GTO observations. In particular, the line fluxes are from 13-19% higher in the ERO observation. 
Thus, both the emission line flux and continuum fiux are lower in the GTO observation nearly five 
months later. 



2.3. Comparison with GHRS and lUE Observations 

In this section, we compare the COS spectra with previous observations taken in the same 
wavelength band. These include a high-resolution GHRS spectrum and two low-resolution lUE 
spectra. While the lUE spectra have a broader wavelength range, similar to the COS spectra, the 
GHRS spectrum covers only the Lya region. We include a comparison of both absorption line 
properties, between the GHRS and COS observations, and the emission line/continuum properties 
between lUE, GHRS, and COS spectra. 



2.3.1. Absorption Line Properties 



To compare the COS ob servations with the G HRS spectrum, we first applied the Galactic 
reddening correction using the ICardelli et al.l (j 19891 ) reddening curve and shifted the spectrum to 
the AGN rest frame. In Figure [H we plot the Lya emission line region for the COS and GHRS 
observations. Several weak absorption components, in addition to Galactic absorption (e.g., from 
S II, Si II, and N V), are detected in both the COS and GHRS observations. 
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The major difference between the COS and GHRS observations of Lya i s that an intrinsic 



absor ber (interpreted as a signature of an outflow emanating from the AGN, e.g.. lBarlow &: Sargent 



19971 . section 1), clearly visible in the GHRS data, is present only at a very low level in the COS 
observations (see Figure|3]). To further quantify the properties of this absorption feature, we include 
the parameters from a Voigt profile fit to the GHRS, individual COS observations, and combined 
COS observations in Table [H We find that the velocity of this feature (~ —4300 km s~^) is roughly 
consistent with the absorber found in the 2000-2001 FUSE observations (~ — 4100kms~^). There 
is no significant difference in the measured FWHM of this feature between the GHRS and COS 
observations, but a slight (~ — lOkms"^) increase in the velocity offset. However, the equivalent 
width of the line decreased by a factor of ~ 5 (from EW ~ 200 mA in 1997 to EW ~ 35 mA m 
2009) between the GHRS and COS observations. We find no difference in this feature between the 
COS ERO and GTO spectra and therefore include in the table measurements from a combined 
COS ERO + GTO fit to this feature (see below for our procedure). 

Additionally, we searched for other absorbers in the GHRS and COS datasets. Since many 
of these lines are quite weak, we maximize our data quality by merging the COS ERO and GTO 
data into a single combined spectrum; the flux in the GTO data was multiplicatively scaled by 
1.4 to match the mean continuum level seen in the ERO data. The flux in both spectra were 
normalized using low-order Legendre polynomials fitted to short regions of the spectr um (~ 10 — 20 



A). F igure [5] compares the absorption lines in the GHRS and combined COS data. IPenton et al 



(j2000l ) identified nine absorption lines, including the absorber described above, (filled stars) in the 
GHRS data at 1222 A < Aobs < 1247.3 A, corresponding to Ly a absorbers at |f | > 1500 kms~^. 



We note three additional absorbers that were not identified in iPenton et al.l (|2000l ) (open stars), 
including one feature which appears in the COS data only. 

We present measurements of the twelve detected Lya absorbers in Table [6l including mea- 
surements of the ~ — 4300 kms~^{cz AGN = —4250 km s~^) listed in Table [5] for a clear com- 
parison with the additional weak absorbers. An adjacent, weaker absorption component to the 
—4300 km s~^ absorber is found at Aobs = 1236.9 A {v = —4100 kms~^). It is likely associated with 
the — 4300 km s~^ feature, as its EW appears similarly diminished in the COS observation. An 
additional weak absorber, coincident with a second 'outflow' component in the FUSE observations , 



also appears to be present in the GHRS and COS spectra at Aobs = 1239.2A. In lOunn et al.l ([2003), 
this second component is shown to arise only in the last of the four FUSE observations, with an 
equivalent width of 0.12A and velocity of — 3680kms~^. The corresponding component in our 
spectra is blended with Galactic N V absorption, but present, with a velocity of ~ — 3550kms~^. 
Additionally, the strength of this feature appears to have decreased between the GHRS and COS 
observations. We estimate an equivalent width for the blended absorption line of 54 it 8mA in 
the 1997 GHRS spectrum and 14 ± 5mA in the 2009 COS spectra. Since this feature is blended 
with N V Galactic absorption, the measurements are uncertain, but it appears that the decrease 
from 1997 to 2009 is consistent with the factor of 5 decrease in the absorbers at — 4250kms~^and 
-4100 kms-^ 
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Of the additional Lya absorbers detected, four absorbers are consistent between GHRS and 
COS (Aobs = I224.2A, 1241. oA, 1249.6A, and 1251.8A). These a bsorbers are hkely interstellar 



features from high- velocity clouds (as described in lPenton et al.ll200Cll ) . An additional four absorbers 



are detected but at a low-significance (< lOo"). Finally, the feature at Aobs = 1243.0, also blended 
with a Galactic N V absorption line, does not appear in the GHRS data but is seen in the combined 
COS spectrum. Thus, this is a fourth weak intrinsic absorber in the AGN spectrum. 

If the weak intrinsic absorption features seen in the FUSE, GHRS, and COS spectra are asso- 
ciated with an ionized outflow, we expect a correlation between changes in the strengths of the lines 
and the luminosity of the source. Particularly, for a factor of five change in the — 4300kms~^ ab- 
sorber, we expect a significant change in the continuum or line flux between the GHRS and COS 
observations. However, this is clearly not the case, as seen in Figure HI Either a change in the 
covering fraction or column density of a cloud in the line of sight is a more likely cause of the 
changes in the intrinsic absorbers. This could be an effect of bulk motions of a cloud/clouds in the 
AGN host galaxy. 



2.3.2. Continuum and Emission-Line Properties 

In order to obtain an estimate of the continuum and Lya line flux in the GHRS observation, 

we fit the Galactic reddening-corrected spectrum with a combination of three Gaussians and a 

power law continuum (oq x (A/1750A)~"). As for the COS spectra, we masked prominent features, 

including the intrinsic absorber. In order to allow some consistency between the COS and GHRS 

spectra, we fixed the power-law index to the average value from the COS spectra (a = 1.30). We 

also constrained the velocity offset and width of the first two Gaussians to the values obtained 

in the COS ERO observation, allowing the amplitude and the parameters of a third Gaussian 

to vary. Our best fit continuum flux density at 1750Awas oq = (0.511 it 0.002) x 10^^'^ erg s^^ 

cm~^ A~^. The flux of the first Gaussian component is 6.76 x 10~^^ ergs s~^ cm~^, while the second 

Gaussian component had a flux of zero. The third Gaussian component had a velocity offset of 
-276.2 ±5.9 km s-\FWHM of 362.6 ±7.1 km s-^ and line flux of 2.57 xl0-i3gj.g g-i ^^-2^ r^^ie 

complex shape of the line, essentially modeled by two Gaussians, was not fit well with this model 
(x^/dof = 3.6), but this simple characterization is sufficient to compare the flux from the GHRS 
spectrum to that in the COS spectra. The total measured line flux is 7.0 x 10~^^ erg s^^ cm~^. 

To compare the GHRS results with the COS observations, we must determine the parameters 
for the continuum and Lya emission using the same narrow band used to fit the GHRS. When we 
re-compute the COS continuum and Lya properties using the same model applied to the GHRS, 
we find ao = (0.554 ± 0.002) x 10"^^ erg s~^ cm"^ A~^ and a total Lya line flux of 8.8 x 10"^^ erg 
s~^ cm~^ for the ERO observation and oq = (0.468 ± 0.003) x 10~^^ erg s~^ cm~^ A~^ and a total 
Lya line flux of 6.4 x 10~^^erg s^^ cm^^ for the GTO observation. These results conflrm what 
a visual inspection of Figure |4] suggests: the ERO observation is the brightest observation in both 
continuum flux and Lya emission. The GHRS measured continuum normalization is 8% lower than 
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the ERO value and 8% higher than the GTO value. The Lya emission from the GHRS observation 
is 20% lower than in the ERO observation and only about 9% higher than in the GTO observation. 
This shows that there is no drastic change in flux between the COS and GHRS observations and 
confirms that the drastic change in the equivalent width of the —4250 km s^^ absorber is probably 
not directly related to the relative luminosity of the source, either in the continuum or the Lya 
emission line. 

Next, we describe our analysis of the lUE spectra and compare their properties to the COS 
and GHRS spectra. The resolution of the lUE spectra {R ~ 200) of Mrk 817 is unfortunately too 
low to search for outflow signatures. It is, however, sufficient to obtain the continuum and emission 
line properties of prominent lines. To determine these parameters, w e first corrected the two lUE 



spectra for the effects of Galactic reddening with the ICardelli et al.l (|l989l ) reddening curve. We 
measured the continuum parameters in the same way as for the original fits to the COS spectra, 
choosing continuum regions free of absorption/emission features with which to fit a power law 
model. We also followed the same procedure to fit the emission lines as for the COS spectra - 
fitting the C IV emission with a broad and intermediate width line and applying the same FWHM 
measurements to the Si IV (+ O IV]), N V, and Lya (two of three Gaussians) emission. The 
measured parameters (power-law parameters and total line flux) are shown in Table [71 In Figure [6l 
we plot the lUE spectra and their best-fit emission line fits for the C IV doublet and Lya -|- N V 
doublet along with the COS ERO spectrum. 

In order to compare our results from the lUE, GHRS, and COS spectra, we plot a light curve 
of the line luminosities and continuum flux densities in Figure [71 Since the GHRS spectra only 
cover Lya, we needed to make a few assumptions to include the GHRS flux for the continuum and 
Lya -|- N V. We assumed that the ratio between the GHRS and COS ERO measurements of the 
continuum and Lya emission from fits of the narrow region from 1185.0-1219.5 A were the same 
as would be measured for the broad-band continuum and Lya -|- N V. Therefore, the GHRS points 
correspond to the COS ERO values from Tables [2l and [3l divided by the factor of 1.08 (continuum) 
and 1.20 (Lya) found from our narrow-band fits. The lUE values are from Table [71 We combine 
Lya and N V to minimize errors associated with separating the blended lines. 

Our results show that the emission line luminosities and continuum flux density levels of the 
lUE observations are lower than those in the GHRS and COS observations. While the Lya -|- 
N V doublet appears to be at a much lower level in the lUE spectra, the lower measurement is 
more likely due to the larger uncertainty associated with the lower-resolution data coupled with 
the complexity of the emission lines and the presence of a strong geo-coronal Lya feature (see 
Figure [6|). Comparing the more reliable measurements of the continuum and C IV doublet, we 
estimate a change of approximately 1.5-1.8 and 1.6-2.3, respectively. Therefore, there appears to 
be no drastic (i.e., order of magnitude) change in the UV continuum or emission lines of Mrk 817 
over the observations - spanning nearly 30 years. 
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3. The X-ray Spectra 

In this section, we describe the results of our analysis of the observations of Mrk 817 from 
XMM-Newton and Swift. We include the observation details and reduction of the X-ray spectra 
in § 13.11 In § 13.21 we include spectral fitting of the XMM-Newton EPIC spectra, with the specific 
goal of determining whether an ionized absorber is present. In § 13.31 we search for variability in the 
spectra, as well as the associated optical/UV photometry, between the XMM-Newton and Swift 
observations. 



3.1. Observations 



Serendipitously, the X-ray observatory XMM-Newton (jJansen et al.ll200ll ) observed Mrk 817 
within 2 weeks of the HST COS GTO observation (the second of the two COS observations). 
Mrk 817 was observed on 2009 December 13 for 14220 s, under the program of PI Brandt (ob- 
servation ID: 0601781401). The data were reduced according to the standard guidelines in the 
XMM-Newton User's Manua0, using the XMM-Newton Science Analysis Software (SAs) version 
9.0.0. We analyzed data from the two types of European Photon Imaging Camera (E PIC) onboard 



XMM -Newton: the two front-illuminated EPIC Metal Oxide Semi-conduc tor (MOS; [Turner et al 



2001 



) CCD arrays and the back-illuminated EPIC pn (|Striider et al.ll200ll ) CCD array The EPIC 
cameras are sensitive in the ~ 0.3 — lOkeV band and have a spectral resolution of E/AE ~ 20-50 
with an angular resolution of ^ 6" FWHM. 

The tasks epchain and emchain were used to reduce the data from the pn and MOS instru- 
ments, respectively. Initially a circular source region of radius 36" was used to extract a source 
spectrum; however, the SAS EPATPLOT task revealed that the observation exhibited "pile-up" in 
the pn and MOSl detectors. Pile-up occurs when two or more X-ray photons deposit charge in a 
single pixel or neighboring pixel during one read-out cycle and are therefore registered as a single 
event with an energy equal to the sum of the photon energies. To correct for pile-up, we followed 
the recommended approach of using an annular source region (inner radius 15", outer radius 45" for 
the pn; the corresponding radii for MOSl were 7.5" and 45"; and the M0S2 data did not display 
pile-up) to extract the spectra for those detectors, until epatplot displayed no strong signature of 
pile-upl£|. Background regions were either chosen to be circles near the source or annuli that exclude 
the central source. Additionally, the background light curves (between 10-12 keV) were inspected 
for flaring, and a comparison of source and background light curves in the same energy range was 
used to determine the portions of the observation in which the background was sufficiently low 
compared to the source. Approximately 75% of the pn observation and 93% of the MOS observa- 



" ^http : //heasarc . nasa . gov/docs/xmni/ sas/USG/ 1 

^See the following thread on dealing with pile-up in the EPIC detectors: 
: //xmm. esac . esa. int/sas/current/documentation/threads/epatplot . shtml] 
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tions were deemed free of flaring; the subsequent spectra were generated from the usable portions 
of each observation. Response matrices and auxiliary files were extracted using the tools rmfgen 
and ARFGEN, and the final spectra were grouped with a minimum of 20 counts per bin using the 
GRPPHA tool. 

The pn spectrum has an exposure time of 7257 s with an average count rate of 4.45ctss~^. 
The MOS spectra have exposure times of 5882 s (MOSl) and 5858s (M0S2), with average count 
rates of 1.95 ctss"! (MOSl) and 3.02 ctss'^ (M0S2). The difference in the MOSl and M0S2 count 
rates is due to the different apertures required since pile-up is detected in MOSl but not M0S2. 
No significant rapid variability is seen in the pn or MOS light curves. 

In addition to the XMM-Newton spectrum, we ir iclude 5 lower signal-to-noise ratio X-ray 
spectra observed with the Swift X-ray telescope (XRT) (IBurrows et al.ll2000l : iHill et al.ll200d ). The 
Swift XRT is sensitive in a similar energy band to the XMM-Newton pn and MOS CCDs (~ 0.3- 
lOkeV). The XRT observations were all downloaded from the High Energy Astrophysics Science 
Archive Research Centefl We use the processed clean event files from the photon counting mode 
to create the spectra. The source spectra were extracted from a circular aperture with a 75" radius, 
centered on the source, using the FTOOL XSELECT. The background spectra were extracted from a 
100" circular region near the position of Mrk 817 and free of additional sources. Standard response 
and ancillary response files from the Swift XRT calibration database were used for the photon 
counting event files with grades 0-12. Details of the XRT observations are included in Table El 



3.2. XMM-Nevirton Spectral Analysis 

Although Mrk 817 is a bright X-ray source, the only broad-band (0.3-10 keV) X-ray spectra 
taken to date are a series of short (< 6ks) snapshots with the Swift XRT (which we discuss in the 
following section) and the 14 ks XMM-Newton observation we now present. The XMM-Newton 
observation of Mrk 817 was obtained as part of a program to stud y the X-ray proper ties of a 
uniform sample of AGN detected in the 22-month Swift BAT catalog (jTueller et al.ll2010l ). In this 
paper, we present a basic analysis of the EPIC spectra - principally to determine whether an X-ray 
outfiow is detected. A more detailed analysis will be presented in Vasudevan et al. (in prep). 



Using XSPEC 12.5 (lArnaudl 11999 ^ . we simultaneously fit the pn, MOSl, and M0S2 X-ray 
spectra with a si mple absorbed powe r law. We accounted for Galactic absorption with the tbabs 
model, using the iDickev k LockmanI (|l99d ^ value of Nh = 1.15 x 10^° cm . We found that the 
addition of intrinsic neutral absorption from the AGN did not improve the fit. The residuals to 
this model did, however, show a soft excess, which we fit with a black body model. This simple 
blackbody -|- power-law (where the power-law model corresponds to F{E) cx E^^^~^)) model is 
statistically a good fit to the data, with reduced /dof = 1.03/399. The best-fit parameters of this 



* ]http : //heasar c . gsf c . nasa . gov/ 1 
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model include F = 2.17 ± 0.06, with a normalization of 6.14 zb 0.19 x 10 ^ photons cm ^ s ^ keV ^ 
at 1 keV, and kT = Q.m^to^Q^i keV, with a normahzation of ^AtW x 10"^ photons cm'^ g-i keV^^ 
at 1 keV. Quoted errors here correspond to the 90% confidence level. The corresponding absorption 
corrected flux in the 0.4-10 keV band is 3.4 x 10~^^ erg s~^ cm~^. 

Sources that exh ibit X-ray signatures of outflows are also known to show UV absorbers (e.g., 
Crenshaw et al.l Il999l ). A simple picture of gas in a single ionization state producing both the 
UV and X-ray absorber features is disproved by simultaneous UV/X-ray grating observations. 
However, these observations have revealed th at the velocity fiel d of the UV and X-ray absorbers 
are similar for a few sources (for instance, see iKaspi et al.l 120021 for details on NGC 3783). Thus, 
the UV absorbers are thought perhaps to arise in higher density knots within the X-ray outflow 
(jKrolik Sz KrisElll995l ). Since no UV absorption signatures were seen in the ionized lines (N V and 
C IV) in the COS GTO spectrum, we expect to find no outflow signatures in the X-ray spectrum 
of Mrk 817. 

X-ray outflow signatures are typically detected in CCD data through the O VII and O VIII 
absorption edges at energies of 0.73 keV and 0.87 keV, respectively. As shown in Figure [H which 
presents the pn spectrum in the soft X-ray band and the ratio of data/model in the same band, there 
is no evidence for the presence of an X-ray warm absorber/outflow. To test this result statistically, 
we added two zedge models, one to account for each of the absorption edges, at the energies of the 
edges. The change in upon adding the edges is very low (Ax^ = 0.58). Calculating errors on the 
optical depth from this model, we find upper limits of rovii < 0.11 and roviii < 0.05 at the 90% 
confidence level, confirming that outflows are not detected in the XMM-Newton X-ray spectrum of 
Mrk 817. 



3.3. Variability 

3.3.1. Long-term variability from the X-ray observations 



While the sensitivity of the Swift XRT observations is too low to search for outflows, the 
XRT spectra are useful for determining changes in both the X-ray flux and overall spectral shape 
over the 2.5 yr period covered with the Swift XRT and XMM-Newton observations. To search for 
va riability between the XMM-Newton and Swift XRT observations, we followed a similar procedure 



to I Winter et al.l (j2008l ). who charact erized X-ray variab ility in 22 Seyferts detected in the Swift 
Burst Alert Telescope AGN survey (jTueller et al.ll2008l ) through simultaneous fits to Swift XRT 
and XMM-Newton spectra. We use the best-fit blackbody + power-law model from the XMM- 
Newton spectra as a baseline model to search for variability between the pn and the five Swift XRT 
spectra. 

In Figure [H we plot the ratio of the data/model of the XMM-Newton pn and five Swift 
XRT spectra, using the best-fit model and parameters described in the previous section for the 
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pn spectrum. Clearly, there is variability in both the spectral shape and the flux between the 
observations. In particular, we find that the two Swift XRT spectra from August 2007 are at a much 
lower flux level than the additional observations. The December 2009 pn spectrum corresponds to 
the highest flux level. 

Allowing for a change in flux between the observations with the addition of a constant model, 
we find a large improvement in the fit with a reduced \^ value of 1.46. Fixing the flux constant 
relative to the pn observation, we find constant values for the Swift XRT observations of 0.50 (2007 
May), 0.60 (2007 July), 0.18 (2007 August 19), 0.19 (2007 August 20), and 0.64 (2009 June). We 
find that allowing an intrinsic neutral absorbing column density to vary does not improve the fit. 

For each of the X-ray spectra, we individually fit a blackbody + power-law model. We fixed 
the temperature of the b lackbody component to 0.10 keV and the Galactic absorption to the 
Dickev &: LockmanI (|l990l ) value. In place of the XSPEC powerlaw model, we used the pegged 
power-law model (pegpwrlw). This model is essentially the same as the power-law model, but the 
normalization is the flux in a defined band (we used 0.3-10 keV). The advantage of this model is 
that it allows us to determine error-bars on the power-law flux. The best-fit parameters for each 
of the X-ray spectra are found in Table [9j The errors on the blackbody normalization are large 
for the Swift observations. However, we find that the XMM-Newton pn spectrum has the highest 
measured blackbody flux normalization. It is unclear whether the higher value in the pn (by a fac- 
tor of ten over the XRT-measured values) is a result of calibration differences or a physical change 
in the emitted spectrum. 

The clearest result of our analysis is that the photon index is strongly correlated with the 
luminosity of the source, as shown in Figure [TOj Using ordinary least-squares linear regression, we 
find that F = (1.13±0.10) log Lo.3_iokeV + (-47.33±4.38), with i?^ = o.94 or Spearman correlation 
coefficient of r = 0.77, corresponding to P = 0.05. Correlations between F and the absorption- 
corrected 2-10 k eV luminosity are perhaps seen in high redshi ft samples (e.g., 



Saez et al.l 120081 ) . but not in low redshift [z < 0.1) samples (e.g., iGeorge et al 



Dai et al 



200C 



2004- 



Winter et al 



20091 ). While low redshift samples do not show the correlation, a correlati on exists through mea- 



surements from mu ltiple observations of individual low redshift sources (e.g.. lMushotzky et al. 



1993 



Winter et al.ll2008l ). The results from our analysis of the Mrk 817 X-ray spectra are consistent with 
previous results, where the photon index becomes steeper for higher luminosities and therefore 
higher accretion rates. 

While broad-band X-ray spectra of Mrk 817 are not available over as long of a tii ne period as 
in th e UV, Mrk 817 was obser ved twice in the soft X-rays (0.1- 2 .4 keV ) with ROSAT (lAschenbach 



1988 



Pfeffermann et al.lll987] ) in the early 1990s, bush et al.l (Il99fil ) find a PSPC count rate of 



0.3992 ± 0.0095 ct s'^ from a spectrum from a 4679 s exposure (Dataset ID: RP701632N00) taken 
1993-11- 06. Additionally, Mrk 817 was included in the ROSAT All-Sky Survey bright source 
catalog (jVoges et al.lll999l '). with a PSPC count rate of 0.102 ± 0.017 cts'^ from a 542 s exposure 
taken 1990-12-06. 
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Using canned ROSAT PS PC response files and an a ssumed 0.1-2.4 keV spectrum of a F 



2.0 power-law absorbed by the iDickev k. LockmanI (|l99d ) value towards Mrk 817, we estimate a 



conversion between the ROSAT count rates and flux (8 x 10 ^^erg s ^ cm ^at lets ^). With this 
conversion, we find that the luminosities in the ROSAT observations are (6.03it0.14) x 10^^ ergss"^ 
(1993 observation) and (1.54 ±0.25) x lO^^ergss"^ (1990 observation). To compare the luminosity 
of Mrk 817 over the nearly 20 yr time span covered from the ROSAT through XMM-Newton 
observations, we calculated the 0.1-2.4 keV luminosity for the Swift XRT and XMM-Newton pn 
observations using the best-fit models in Table El 

In Figure [TT| we plot a light curve of the observed X-ray luminosity in the 0.1-2.4 keV band 
constructed from the ROSAT, Swift XRT, and XMM-Newton observations of Mrk 817. Variability 
is clear, with the 1990 ROSAT observation being 39 times less luminous than the 2009 XMM- 
Newton observation. A factor of four variability is seen between the two ROSAT observations, 
which are separated by three years. Additionally, in the Swift XRT observations, which also span 
approximately three years, there is a factor of four variability. Thus, the X-ray spectra of Mrk 817, 
which span ~ 20 years, show a higher level of variability than was observed in the UV spectra (from 
lUE, GHRS, and COS). 



3.3.2. Variability from the simultaneous Optical/UV/X-ray data 
To test whether the variability in the X-ray band is correlated with opt ical/UV variability in 



Mrk 817, we utilized XMM-Newton Optical Monitor fOM: lMason et al.ll200ll ) and Swift Ultraviolet 



Optical Telescope (uvotIiI) data that were taken simultaneously with the X-ray spectra. To extract 
the OM fluxes, we ran the standard SAS omichain script to produce processed OM images. The 
fluxes were obtained from the omdetect algorithm, which extracts the flux for point-like sources in 
the images. OM images were taken in both the UVWl (Ac = 2910 A) and UVM2 (Ac = 2310 A) 
filters for our XMM-Newton observation. UVOT observations were available for three of the five 
Swift observations of Mrk 817 (00036620001, 00036620002, and 00037592001). To determine the 
flux of our source in the Swift UVOT observations, we ran the Swift task uvotsource on the 
processed UVOT images, downloaded from HE AS ARC. We used a source region centered on Mrk 
817 with a 5" circular aperture. A background region was selected on the image in an area free 
from bright sources. In Table [T0| we list the measured fluxes from each of the available filters. 
Since Mrk 817 is a nearby face-on spiral, the extracted OM/UVOT fluxes are a combination of 
both the AGN and stellar bulge emission. 

In Figure m we plot the observed fluxes obtained from UVOT and OM U, UVWl, UVM2, 
and UVW2 filters. We also plot the continuum values measured from our COS observations and 
extended to 4000 A. We find that the UVOT and OM fluxes are well-matched to the slope and flux 



^See the UVOT website at: [http: //www. swift .psu. edu/uvot/| 
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of the COS observations. In particular, the OM values, taken within two weeks of the COS GTO 
observation, are < 10% higher than the estimated continuum from COS. As shown in the figure, 
the scale of optical/UV variability over the Swift, COS, and XMM-Newton observations is minimal 
(~ 30%) over the 2.5 year time frame. Comparing the Swift UVOT U-band fluxes (the only filter 
for which we have data for the three Swift observations), we find that the highest flux observation 
is (33 lb 6)% higher than the lowest flux observation. The X-ray variability is comparable, with 
a (27 ± 5)% change between the highest and lowest flux observation. However, the observation 
with the highest X-ray flux does not correspond to the observation with the highest U-band flux. 
Therefore, there is no obvious correlation between the X-ray and U-band flux. 

The availability of simultaneous, multi-band data allows the construction of multi-epoch SEDs 
for Mrk 817 along with simple model fits to the data. This illustrates how the value of aox, which is 
the optical-to-X-ra-y spec tral slope, or X-ray loudness, defined as aox = —0.384 log (F2kev/-^ :) 

Zo(JuJ\ 



(ITananbaum et al.lll979l'). varies between different observations. We follow the approach outlined 



m iVasudevan et al.l (120091 ) for generating the SEDs, using the ftool flx2xsp to convert the 
XMM-Newton OM and Swift UVOT photometry into data files suitable for fitting in xspec. 
We constructed SEDs for the three Swift observations with IDs 00036620001, 00036620002 and 
00037592001 and the new XMM-Newton observation. A simple disk plus absorbed power-law 
model was then fit to the optical/UV through X-ray data points (full model used: diskpn + 



wabs(zwabs(bknpower))), with the model parameters constrained as specified in lVasudevan et al 



20091 ). The no rmalization o f the diskpn component was constrained using the reverberation mass 



estimate from iBentz et al.l ( 20091 ) and Galactic neutral hydrogen was included via the wabs com- 
ponent. 

The full SEDs for the four observations under discussion, along with their respective model fits 
(with any absorption removed to show the intrinsic SED model), are shown in Figure[T3l The model 
fits were then used to determine the value of aox, by computing the monochromatic luminosities 
at rest-frame 2500 A and 2keV. We obtain the following values of aox for the observations: XMM- 
Newton: -1.353+[!:JJ'i'?, Swift 00037592001: -1.51+|]:?^ and Swift 00036620001: -1.5ltofT The 
computed aox value for the observation Swift 00036620002 is omitted due to its high uncertainty, 
since UVOT data was available in only one filter (U). The XMM-Newton observation has the highest 
X-ray luminosity (5.9 x 10'*^ ergs s~^ in the 0.3-10 keV band) and the smallest value of aox ■ The two 
Swift observations with multiple UVOT filter observations are consistent with having the same aox 
values. The relationship between aox-^x is not entirely clear, however, since the Swift observations 
have the same aox value but vary in luminosity by 26.9 it 4.6%. Additionally, the XMM-Newton 
observation is only 18.9 it 4.3% more luminous than the Swift 00037592001 observation, but has a 
significantly different aox value. 
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4. Conclusions 

In this paper, we present an analysis of several UV and X-ray spectra of the Seyfert 1.5 Mrk 
817. The UV spectral analysis includes new high-resolution spectra from the Cosmic Origins Spec- 
trograph, as well as archived spectra from lUE and HST GHRS. These observations span nearly 
30 years and allow us to search for spectral variability on decade time scales. The X-ray spectral 
analysis includes a spectrum from XMM-Newton, taken within two weeks of the second COS obser- 
vation, as well as five lower signal-to-noise ratio spectra with Swift XRT. These observations allow 
us to probe X-ray spectral variability on the 2.5 year time scale covered. In addition, luminosities 
from ROS AT allow us to determine variability in the soft X-ray band over a 20 year period. Finally, 
an analysis of the OM/UVOT observations, which were simultaneous with the X-ray observations, 
allows us to determine variability in the UV/X-ray SED. In this section, we summarize and discuss 
the results of our analysis. 



4.1. Spectral Shape and Fundamental Properties 



The HST COS observations of Mrk 817 exhibit the broad UV emission lines (Lya, N V, Si IV 
-|- O IV], and C IV) characteristic of an optical broad-line Seyfert. Our analysis of the ERO and 
GTO observations, separated by ~ 5 months, finds similar power-law slopes {Fx oc A"^'^'') and 
FWHMs for the emission lines (e.g., we find a broad component of ~ 3000 km s~^ and a narrower 
component of ~ 1000 km s~^ for the C IV emission). The C IV emission is blue-shifted by values 
similar to other AGNs (blue-shifted by 50-340 kms~^), but there is no correlation between the 
shifts in the C IV/Si IV/N V/Lya emission components. Further, these velocity shifts change 
between observations. 

No strong intrinsic absorption lines (i.e., Lya, N V, C IV) are detected in the COS observations. 
However, four high-velocity cloud absorbers and four intrinsic Lya absorption lines are detected 
in the GHRS and COS spectra (as well as four absorbers detected at low-significance, see § I2.3.ip . 
Corresponding features are not detected in the more ionized lines. 

The X-ray spectra of Mrk 817 are well-characterized by a blackbody + power-law model. The 
temperature of the blackb ody component, us e d to characterize th e soft excess, of kT 0.1 keV, is 



typic al of Seyfert Is (e.g., iMateos et al.ll2005l : IWinter et al.ll2009l ) and quasars (ICierlihski fc Done 
20041 ). However, the power-law slope is steeper than the average photon index of similar AGN. We 
find r > 2.0 in all but the lowest flux observations; while typical values for local Seyfert Is detected 
in the hard X-rays with Swi ft's BAT (of which , Mrk 817 is included in the 22-month Swift BAT 
catalog) are F = 1.78 ± 0.24 (jWinter et al.ll2009l ). Other than the steep slope, there are no unusual 
features in the spectrum. Additionally, we find no evidence for either neutral or ionized absorption. 



Using the flux and continuum properties derived from our spectral analysis, we determine 
fundamental properties of Mrk 817 to determine how it compares to local Seyferts. The X-ray 
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2-10 keV luminosity, which we compute using the 2-10 keV flux of 1.45 x 10~^^ erg s~^ cm~^ from 
the XMM-Newton spectrum, can be used to estimate the bolo metric luminosity o f the AGN. Using 
a bolometric correction of 35, typical of unabsorbed AGN (jBarger et alj |2005| ) and within the 
range of parameters determined from our SED fits to the simultaneous optical/UV/X-ray data 
from Swift and XMM-Newton, we estimate Lf,oi = 9.58 x 10^^ ergs~^. This value is consistent with 
local {{z) = 0.03) broad-line Seyferts selected in the 14-195 keV band with the Swift BAT detector, 
where the average bolometric lumino sity is 7.4 x 10^^ ergs~^, with a range of values from ~ 9 x 10^^ 
- 6 X 10^^ ergs"^ (jwinter et a llboiol) . The black hole m ass of Mrk 817, log M/M© = 7.69 ±0.07 as 
deriv ed from reverberation mapping (jBentz et al.ll2009l ). is also typical of loca l AGN (Woo &: Urrv 



2OO2I ) and the BAT AGN sample in particular ( (log M/Mg 



Winter et al 



7.87 ± 0.66; 

However, the estimated accretion rate {Li,oi/ LeiMi where Lsdd = {M/Mq) x 1.3 x 10"^*ergs~^) of 
0.14 is higher than the average accretion rate of 0.034 (with a range of va lues from ~ 0.001-0.5) 
determined for broad-line Seyferts in the BAT sample (jWinter et al.ll2010l ). Thus, Mrk 817 is in 
many ways (mass, luminosity, UV spectral shape) a typical broad-line Seyfert. Since the accretion 
rate is believed to determine the har d X-ray spectral slop e in AGNs, with high spectral slopes 



201 



3). 



correlating with high accretion rates (jShemmer et al.ll2006l ). it is not surprising that Mrk 817 has 



both a slightly higher photon index and accretion rate than the average local Seyfert 1. 



4.2. Intrinsic Absorption 

Four intrinsic absorption features were detected embedded in the broad Lya emission of Mrk 
817. These features include a —4250 km s~^, — 4100kms~^, — 3550 km s~^, and —2600 km s~^ absorber. 
The strength of these features vary between the 1997 GHRS and 2009 COS observations, with the 
ratio of EW between GHRS/COS being 5.7 ± 1.3, 2.2 ± 1.1, 3.9 ± 3.0, and ~ 0.3, respectively. 
The width of the absorption features are consistent between the GHRS and COS observations. 
No intrinsic C IV AA1548. 20,1550. 77 and N V AA1238.82,1242.80 absorption features are detected 
at the same velocities. However, given the weakness of the intrinsic Lya absorbers in the COS 
spectra (EW ~ 35 mA), this absence is expected. For similar ionic column densities, the ratio 
of optical depth between two lines should be proportional to the ratio of their /A values (where 
/ is the oscillator strength and A is the wavelength of the transition). Under these conditions, 
r(Lya)/r(C IV A1548)=1.7 and r(Lya)/r(N V A1238)=2.6. For a photo-ionized plasma with a 
fixed total hydrogen column density, the column densit y of each ion is a strong function of the 



ionization parameter U. However, as shown in figure 2 of lArav et al.l (j2007l ). for solar abundances, 
the column densities of C IV and N V are always smaller than that of H I, irrespective of the 
value of U. Combined with the stronger expected optical depth of the Lya line, a weak absorption 
feature of the latter should not show related features in C IV and N V. The same figure shows that 
under similar conditions O VI absorption is expected for log U ~ —0.5, where for solar abundances 
the O VI column density is 8 times larger than that of H I, while r(Lya)/r(0 VI A1032)=3.7 for 
the same column densities. Under these conditions, we expect r(0 VI A1032) to be twice that of 
Lya. This may explain why O VI and Ly/3 features are detected in the same FUSE data (For the 
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same column densities, r(Ly/3)/r(0 VI A1032)=0.6.). The photo-ionization curves also show that 
r(0 VI A1032)/r(Lya) drops below 1 for log [/ > 0.4 and logC/ < —1. Similarly, with such a weak 
Lya absorption feature, we do not expect to detect the O VII and O VIII absorption edges in the 
near-simultaneous X-ray data. 

The strongest GHRS Lya absorber at —4250 kms^^, which is 5 times weaker in the COS 
spectra, is likely associated with th e strongest O VI absorber seen in the four FUSE observations 
from 2000-2001 (jPunn et ahlbood l. As we found for the GHRS/COS observations, the O VI 



absorber, whose velocity changes slightly from —4198 km s ^ to —4144 km s ^ throughout the ob- 
servations, also varies in EW (from 0.60 A to 0.33 A). Both the strength and velocity of the O VI 
absorber decreased over time - however, this decrease was not correlated with changes in the Far- 
UV luminosity of Mrk 817. Additional weak O VI absorption lines were also found to appear 
throughout the set of four FUSE observations, with no obvious correlation to the UV luminosity. 
In a similar way, a comparison of the COS observations with the GHRS observation shows no 
correlation between UV luminosity (in the continuum or emission lines) and the equivalent width 
of the absorber. 

Mrk 817 is not alone in exhibiting variable absorbers. Va riability in UV absorption lines has 
beeii detected in several S eyfert Is, includi ng NGC 4151 



1985 
1993 



Kraemer et al 



Crenshaw et al.l 



20091 ^. and NGC 3783 (e.g.. 



Ulrich & Boisson 


1 


ad), NGC 5548 (e.g.. 


S 


hull & Sachs 


(e.g., 


Mar an et al. 


1996; 


Gabel et al. 


2005 


) . The cause 



of variability of intrinsic absorbers is attributed to a change in ionization state or shielding of the 
photo- ionizing continuum of an outflowing wind or bulk motions of an absorbing cloud/clouds. 

A potential cause of the change in th e absorbers ii i Mrk 817 is from bulk motions of clouds 
along our line-of-sight. As pointed out in iDunn et al.l (j2008l ). weak absorption - as in the O VI 
absorbers in Mrk 817 - is seen in spectra where this is the case. Further, the fact that there is little 
change in the UV luminosity between the GHRS and COS spectra, while the strength of the Lya 
absorber decreases significantly, also favors this model. This would suggest that ionization is not 
causing a change in the absorber. Based on the assumption of an absorbing cloud moving radially 
out of the line of sight in the 12.5 years between the GHRS and COS observations, this implies 
that the cloud moved a distance about 10^ ''cm (<C Ipc). 

Based on the X-ray measured upper limits on the O VII and O VIII absorption edges, we can 
estimate an upper limit on the column density of the absorbing cloud to compare w i th UV -derived 
measurements. Using the cross-sections {a) for O VII and O VIII from lverner et all ( 1996 ) and the 
relation that r = (xNo.ioni we find upper limits on the column density as Nqvii ^ 4.0 x 10^'' cm^^ 
and Noviii ^ 4.7 x lO^^cm"^. The typical range of column densiti es in C IV is well bel ow the 



limit we find for O VII and O VIII, with Nciv = (0.1-14) x 10^^ cm'^ (jCrenshaw et al.lll999l ). This 
X-ray derived limit is also much higher than the estimated H I column from the Lya absorption 
line (~ 4.3 x 10^^ cm~^ in the GHRS and ~ 7.1 x 10^^ cm~^ in COS, assuming a fully covering 
absorber). Therefore, it is possible that an outflow is present, but it has a column density too low 
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to be detected in the CCD spectra analyzed. 



4.3. Variability 

With available observations of Mrk 817 from multiple observatories spanning decades, we were 
able to probe variability in flux and spectral shape across the UV and X-ray bands. In the UV, from 
analysis of the lUE, GHRS, and COS spectra, we do not find drastic changes in either the line or 
continuum flux of Mrk 817. As with other AGN, the c ontinuurn flux is correlated with t he emission 
line flux of C IV in accordance with the Baldwin effect ( Baldwinlll977 : Osmer et al. 1994 ). Measured 



changes in UV flux are limited to a factor of < 2.3 between the minimum and maximum measured 
values over the 30 year time-span probed. Additionally, in support of our spectral analysis, we find 
that the UV photometry from XMM-Newton OM and Swift UVOT observations are at a similar 
level as the UV continuum measured from the COS spectra. 

Variability in the X-ray spectra is much more pronounced than in the UV. In particular, we find 
that the X-ray spectrum of Mrk 817 varies in both luminosity and spectral shape. The photon index 
measured from a power-law component varies from T ~ 1.5-2.1 between the five Swift XRT and 
the XMM-Newton spectra, which span 2.5 years. The variability in the spectral slope is correlated 
strongly with the observed X-ray luminosity, where the steepest slopes are observed at the highest 
X-ray luminosities. 

Two ROSAT flux measurements from the early 1990s have allowed us to search for long- 
term variability in the soft X-ray band luminosities from the ROSAT, Swift, and XMM-Newton 
observations. We find that the X-ray luminosity of Mrk 817 changes markedly. Shorter term 
variability is observed in the Swift XRT observations, with a factor of four change in variability 
over ~ 2 years. On the long term, we find that the 1990 ROSAT luminosity is ~ 40 times smaller 
than the luminosity measured from the 2009 XMM-Newton observation. The factor of 40 variability 
seen over 20 years in the X-ray observations contrasts with the factor of 2.3 seen in the UV spectra. 
Thus, Mrk 817 is clearly more variable in the X-ray than UV band. 

To search for a possible relationship between the UV and X-ray luminosities, we analyzed 
simultaneous optical/UV/X-ray data from XMM-Newton and Swift. We found that the optical/UV 
fluxes were at similar levels in the four observations with OM/UVOT photometry. The optical- X- 
ray SED of Mrk 817 changes between observations, as indicated by aox- However, these changes 
are dominated by the change in X-ray luminosity /spectral shape. Therefore, it appears that the 
UV is not correlated with the X-ray and that even in our simultaneous UV/X-ray data, we find 
more variability in the X-ray than the UV. Unfortunately, however, there is no UVOT data for 
the two Swift observations with the lowest X-ray luminosities, making it impossible to determine 
whether the UV luminosity is low when the X-ray luminosity is at a much lower level than sampled 
with the available simultaneous data. 



Finally, our analysis of the UV and X-ray observations suggests that more data are needed 
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before we can rule out the possibility that the Lya absorbers are the result of an ionized outflow. 
While wc did not find an X-ray O VII or O VIII absorption edge strongly detected in our data, 
it is possible that a high signal-to-noise grating observation could reveal ionized X-ray absorp- 
tion/emission at a lower column density than we could probe with the CCD spectra. Additionally, 
while the UV luminosity does not change significantly, it is possible that the X-ray spectrum during 
the 1997 GHRS observation was significantly different than in 2009, when the COS spectrum was 
taken. Further, our results of both a lack of a correlation between the UV and X-ray luminosi- 
ties and a nearly constant UV luminosity over 30 years while the X-ray luminosity changes by an 
order of magnitude, coupled with the possibility that changes in ionization of an outflow are the 
result of changes in the X-ray emission, suggest that the factor of 5 change in the Lya absorption 
is potentially explained by a change in the X-ray spectrum. Therefore, future multi-wavelength 
observations, with COS and X-ray grating data from Chandra/XMM-Newton, of Mrk 817 are vital 
to inivcil both the true nature of the absorbers and any potential connection between the UV and 
X-ray emission. 

The authors acknowledge support through NASA grant HST-NNX08AC14G, in support of 
the COS GTO program (PI Green). LMW acknowledges support through NASA grant HST-HF- 
51263. 01-A, through a Hubble Fellowship from the Space Telescope Science Institute, which is 
operated by the Association of Universities for Research in Astronomy, Incorporated, under NASA 
contract NAS5-26555. WNB and RVV acknowledge support from NASA grant NNX09AQ02G and 
NASA ADP grant NNX10AC99G. Some of the data presented in this paper were obtained from the 
Multimission Archive at the Space Telescope Science Institute (MAST). STScI is operated by the 
Association of Universities for Research in Astronomy, Inc., under NASA contract NAS5-26555. 
Support for MAST for non-HST data is provided by the NASA Office of Space Science via grant 
NNX09AF08G and by other grants and contracts. This work utilizes observations obtained with 
XMM-Newton, an ESA science mission with instruments and contributions directly funded by ESA 
Member States and NASA. 

Facilities: HST(), Swift(), XMM(), IUE() 

REFERENCES 

Arav, N. et al. 2007, ApJ, 658, 829 

Arnaud, K. 1996, Astronomical Data Analysis Software and Systems V, 101, 5 
Aschenbach, B. 1988, Appl. Opt., 27, 1404 
Baldwin, J. A. 1977, ApJ, 214, 679 

Barger, A. J., Cowie, L. L., Mushotzky, R. F., Yang, Y., Wang, W.-H., Steffen, A. T., & Capak, 
P. 2005, AJ, 129, 578 



- 22 - 



Baxlow, T. A., & Sargent, W. L. W. 1997, AJ, 113, 136 
Baskin, A., k Laor, A. 2005, MNRAS, 356, 1029 

Bentz, M. C, Peterson, B. M., Netzer, H., Pogge, R. W., & Vestergaard, M. 2009, ApJ, 697, 160 
Bromage, G. E. et al. 1985, MNRAS, 215, 1 

Burrows, D. N. et al. 2000, in Presented at the Society of Photo-Optical Instrumentation Engineers 
(SPIE) Conference, Vol. 4140, Society of Photo-Optical Instrumentation Engineers (SPIE) 
Conference Series, ed. K. A. Flanagan &; O. H. Siegmund, 64-75 

CardeUi, J. A., Clayton, G. C, & Mathis, J. S. 1989, ApJ, 345, 245 

Collins, J. A., ShuU, J. M., & Giroux, M. L. 2003, ApJ, 585, 336 

Constantin, A., & Shields, J. C. 2003, PASP, 115, 592 

Crenshaw, D. M., Kraemer, S. B., Boggess, A., Maran, S. P., Mushotzky, R. F., & Wu, C.-C. 1999, 
ApJ, 516, 750 

Crenshaw, D. M., Kraemer, S. B., Schmitt, H. R., Kaastra, J. S., Arav, N., Gabel, J. R., & Korista, 
K. T. 2009, ApJ, 698, 281 

Dai, X., Chartas, G., Eracleous, M., & Garmire, G. P. 2004, ApJ, 605, 45 

Danforth, C. W., Keeney, B. A., Stocke, J. T., ShuU, J. M., & Yao, Y. 2010, ArXiv e-prints 

Dickey, J. M., & Lockman, F. J. 1990, ARA&A, 28, 215 

Dixon, W. et al. 2010, Cosmic Origins Spectrograph Instrument Handbook, 2nd edn., STScI, 
Baltimore 

Dunn, J. P., Crenshaw, D. M., Kraemer, S. B., & Gabel, J. R. 2007, AJ, 134, 1061 
Dunn, J. P., Crenshaw, D. M., Kraemer, S. B., & Trippe, M. L. 2008, AJ, 136, 1201 
Falco, E. E. et al. 1999, PASP, 111, 438 
Gabel, J. R. et al. 2005, ApJ, 631, 741 

George, I. M., Turner, T. J., Yaqoob, T., Netzer, H., Laor, A., Mushotzky, R. F., Nandra, K., & 
Takahashi, T. 2000, ApJ, 531, 52 

Gibson, B. K., Giroux, M. L., Penton, S. V., Stocke, J. T., ShuU, J. M., & Tumlinson, J. 2001, AJ, 
122, 3280 

Gierlinski, M., & Done, C. 2004, MNRAS, 349, L7 



-23- 



Green, J. C, Science, C, & Instrument Development Teams. 2010, in Bulletin of the American 
Astronomical Society, Vol. 41, Bulletin of the American Astronomical Society, 501 — h 

Hill, J. E., Zugger, M. E., Shoemaker, J., Witherite, M. E., Koch, T. S., Chou, L. L., Case, T., & 
Burrows, D. N. 2000, in Presented at the Society of Photo-Optical Instrumentation Engineers 
(SPIE) Conference, Vol. 4140, Society of Photo-Optical Instrumentation Engineers (SPIE) 
Conference Series, ed. K. A. Flanagan & O. H. Siegmund, 87-98 

Ilic, D., Popovic, L. C., Bon, E., Mediavilla, E. G., & Chavushyan, V. H. 2006, MNRAS, 371, 1610 

Jansen, F. et al. 2001, A&A, 365, LI 

Kaspi, S. et al. 2002, ApJ, 574, 643 

Kraemer, S. B. et al. 2006, ApJS, 167, 161 

Krolik, J. H., & Kriss, G. A. 1995, ApJ, 447, 512 

Maran, S. P., Crenshaw, D. M., Mushotzky, R. F., Rcichcrt, G. A., Carpenter, K. G., Smith, A. M., 
Hutchings, J. B., & Wcymami, R. J. 1996, ApJ, 465, 733 

Mason, K. O. et al. 2001, A&A, 365, L36 

Mateos, S., Barcons, X., Carrera, F. J., Ceballos, M. T., Hasinger, G., Lehmann, I., Fabian, A. C, 
& Streblyanska, A. 2005, A&A, 444, 79 

Mushotzky, R. F., Done, C, & Pounds, K. A. 1993, ARA&A, 31, 717 

O'DonneU, J. E. 1994, ApJ, 437, 262 

Osmer, P. S., Porter, A. C, & Green, R. F. 1994, ApJ, 436, 678 

Osterman, S. et al. 2010, in Bulletin of the American Astronomical Society, Vol. 41, Bulletin of the 
American Astronomical Society, 502 — h 

Penston, M. V. et al. 1981, MNRAS, 196, 857 

Penton, S. V., Stocke, J. T., & Shull, J. M. 2000, ApJS, 130, 121 

Peterson, B. M. et al. 2004, ApJ, 613, 682 

Pfeffermann, E. et al. 1987, in Society of Photo-Optical Instrumentation Engineers (SPIE) Confer- 
ence Series, Vol. 733, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference 
Series, 519—1- 

Predehl, P., & Schmitt, J. H. M. M. 1995, A&A, 293, 889 

Rush, B., Malkan, M. A., Fink, H. H., k Voges, W. 1996, ApJ, 471, 190 



-24- 



Saez, C, Chartas, G., Brandt, W. N., Lehmer, B. D., Bauer, F. E., Dai, X., k Garmire, G. P. 
2008, AJ, 135, 1505 

Scott, J. E., Kriss, G. A., Brotherton, M., Green, R. F., Hutchings, J., Shull, J. M., & Zheng, W. 
2004, ApJ, 615, 135 

Shang, Z., Wills, B. J., Wills, D., & Brotherton, M. S. 2007, AJ, 134, 294 

Shemmer, O., Brandt, W. N., Netzer, H., Maiolino, R., k Kaspi, S. 2006, ApJ, 646, L29 

Shull, J. M., & Sachs, E. R. 1993, ApJ, 416, 536 

Striider, L. et al. 2001, A&A, 365, L18 

Tananbaum, H. et al. 1979, ApJ, 234, L9 

Telfer, R. C., Zheng, W., Kriss, G. A., & Davidsen, A. F. 2002, ApJ, 565, 773 
Tueller, J. et al. 2010, ApJS, 186, 378 

Tueller, J., Mushotzky, R. P., Bartliclmy, S., Cannizzo, J. K., Gehrels, N., Markwardt, C. B., 
Skinner, G. K., & Winter, L. M. 2008, ApJ, 681, 113 

Turner, M. J. L. et al. 2001, A&A, 365, L27 

Ulrich, M. H., & Boisson, C. 1983, ApJ, 267, 515 

Vasudevan, R. V., Mushotzky, R. F., Winter, L. M., & Fabian, A. C. 2009, MNRAS, 399, 1553 
Verner, D. A., Ferland, G. J., Korista, K. T., & Yakovlev, D. G. 1996, ApJ, 465, 487 
Voges, W. et al. 1999, A&A, 349, 389 
Wilkes, B. J. 1986, MNRAS, 218, 331 

Winter, L. M., Lewis, K. T., Koss, M., Veilleux, S., Keeney, B., & Mushotzky, R. F. 2010, ApJ, 
710, 503 

Winter, L. M., Mushotzky, R. F., Reynolds, C. S., & Tueller, J. 2009, ApJ, 690, 1322 
Winter, L. M., Mushotzky, R. F., Tueller, J., & Markwardt, C. 2008, ApJ, 674, 686 
Woo, J.-H., k Urry, C. M. 2002, ApJ, 579, 530 



This preprint was prepared with the AAS IM^^X macros v5.2. 



-25- 



Table 1. Summary of Ultraviolet Observations 



Observation 


Date 


Exp Time (ks) 


Spectral Coverage 


A/AAt 


COS ERO 


2009-08-04 


3.4 


1135-1795 A 


18000 


COS GTO 


2009-12-28 


3.0 


1135-1795 A 


18000 


GHRS 


1997-01-12 


25.0 


1223-1259 A 


18000 


lUE SWP 15447 


1981-11-07 


16.2 


1120-2010 A 


200 


lUE SWP 17442 


1982-07-18 


16.2 


1120-2010 A 


200 



^Approximate spectral resolution at 1200 A. 



Table 2. Power law Best-fit Parameters 



Parameter 


COS ERO 


COS GTO 


ao 


0.63 ± 1.02 


0.45 ± 1.02 


a 


1.30 ±0.05 


1.22 ±0.07 



Parameters of a simple power law fit, of the form a,o x (A/1750 A) where qq is the continuum 
flux density at 1750 A (in units of 10~^^ erg s~^ cm~^ A~^) and a is the power law index. 
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Table 3. Emission Line Best-fit Parameters for the COS ERO spectrum 



Component 


v(off) (kms-l) 


FWHM (kms-i) 


EW (A) 


Log Flux 






CIV 1550.774 A 










;^2/dof = 9402.9/9491 






1 


-337.6 ± 15.1 


3145 ±31 


15.97 ±0.14 


-11.93 ± 0.01 


2 


-126.1 ± 15.4 


1051 ± 30 


4.04 ± 0.17 


-12.53 ± 0.02 






Si IV 1402.77 A 










X^/do{ = 21104.7/11489 






1 


381.3 ±22.6 




4.45 ± 0.04 


-12.43 ± 0.01 


2 


432.2 ±55.5 




0.45 ± 0.02 


-13.15 ±0.01 






NV 1242. 80A 










xVdof = 20224.6/7807 






1 


27.2 ±8.6 




7.69 ± 0.04 


-12.12 ± 0.01 


2 


-100.0 ±2.9 




1.32 ±0.02 


-12.89 ± 0.01 






Lyo 1215. 67A 






1 


-905.0 ±9.8 




48.52 ±0.47 


-11.31 ± 0.01 


2 


800.0 ± 10.0 




12.82 ±0.17 


-11.89 ±0.01 


3 


-556.4 ±22.0 


1900 ± 26 


17.20 ±0.38 


-11.76 ±0.01 



The results of multiple Gaussian fits, to the indicated UV emission lines, are shown. The FWHM for Si IV, N V, and the first 
two components of hya are fixed to the values from the CIV fits. The logarithm of the integrated line flux is recorded for 

each component (with units of ergs cm"^). We fit the doublets (C IV, Si IV, and N V) by fixing the width between the two 
lines of the doublet to the known separation of the lines. The flux and FWHM of each component in the doublet were tied 

together, with the flux of the blue line in the doublet fixed to be twice that of the red line (the ratio of the oscillator strengths 

of blue/red is ^ 2 : 1 for each of these doublets). The total flux for each doublet, then, is three times the total flux for the red 

line listed in the Table. 



Table 4. Emission Line Best-fit Parameters for the COS GTO spectrum 



Component 


v(ofT) (kms 


) FWHM (kms-i) 


EW (A) 


Log Flux 






CIV 1550.774 A 










x2/dof = 8972.2/9491 






1 


-53.5 ±17.5 


3481 ± 33 


18.05 ±0.14 


-12.03 ±0.01 


2 


-270.0 ± 12.4 


1030 ± 19 


5.00 ±0.14 


-12.59 ±0.01 






Si IV 1402.77 A 










;^2/dof = 14262.0/11491 






1 


765.8 ±27.5 




5.38 ± 0.05 


-12.50 ±0.01 


2 


232.2 ±48.0 




0.65 ± 0.02 


-13.41 ±0.01 






NV 1242.80 A 










xVdof = 11948.1/7807 






1 


-216.1 ± 14.6 




9.67 ±0.08 


-12.18 ±0.01 


2 


-100.0 ±4.9 




1.81 ± 0.04 


-12.91 ±0.01 






Lyo 1215.67A 






1 


-1154.2 ±8.8 




57.57 ±0.30 


-11.40 ±0.01 


2 


800 ± 5.8 




18.07 ±0.37 


-11.90 ±0.01 


3 


-56.7 ±28.1 


1833 ± 34 


16.90 ±0.39 


-11.93 ±0.01 



The notes to this table are the same as for Table |3] 
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Table 5. Intrinsic Lya Absorber Measurements 



Observation 


v(sys) (kms ^) 


FWHM (kms-i) 


EW (mA) 


logNm (cm 2) 


GHRS 


-4269 ± 3 


125 ±8 


199 ± 10 


13.63 ±0.02 


COS ERO 


-4332 ± 20 


113 ± 17 


39 ±10 


12.9 ±0.1 


COS GTO 


-4322 ± 20 


< 133 


53 ±10 


13.0 ±0.1 


Combined COS 


-4288 ± 10 


110 ± 17 


37 ±5 


12.85 ±0.08 



The parameters of Voigt profile fits to the strongest absorption feature in the GHRS and COS 
spectra include the systemic velocity of the absorber, the full-width at half maximum, equivalent 
width, and H I column density (assuming a fully covering absorber) . 



Table 6. Intrinsic Lya Absorbers 



^obs 


czagn'' 


GHRS EW 


GHRS FWHM 


COS EW 


COS FWHM 


Absorber 


(A) 


( km s^'-) 


(mA) 


( km s^'-) 


(mA) 


(kms-i) 


Notes 


1223.5'" 


-7410 


19 ±10 


43 ±23 


< 6 




low-significancet 


1224.2^' 


-7250 


130 ± 10 


63 ±13 


150 ±6 


75 ±8 


consistent 


1234.7'' 


-4660 


33 ±15 


83 ±17 


< 5 




low-significancet 


1236.3'' 


-4250 


200 ± 10 


125 ±8 


35 ±5 


110 ±8 


variable 


1236.9'' 


-4100 


28 ±5 


40 ±8 


13 ±3 


55 ± 20 


weaker component 


1239.2'' 


-3550 


54 ±8 


122 ± 17 


14 ±5 


60 ± 13 


Galactic N V blend 


1241.0'' 


-3090 


34 ±5 


40 ±10 


34 ±3 


47 ± 10 


consistent 


1243.0 


-2600 


< 5 




20 : 


58 : 


Galactic N V blend 


1245.4'' 


-2010 


18 : 


75 : 


< 5 




low significancet 


1247.3'' 


-1540 


30 ±6 


107 ± 17 


< 5 




low significancet 


1249.6 


-980 


66 ±5 


85 ±8 


64 ±5 


73 ±8 


consistent 


1251.8 


-440 


14 ±4 


25 ±8 


8±5 


25 : 


weak but consistent 



^'Offset velocity from assumed AGN velocity of 9341 kms ^. 

''Absorber listed in Penton et al. (2000). 

t Low significance is defined as significance < 10a. 
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Tablc 7. Best-fit Parameters for tfie lUE Spectra 



Parameter 


SWP 15447 


SWP 17442 




0.35 ± 1.01 


0.41 ± 1.01 




0.55 ±0.07 


0.98 ± 0.05 


C IV doublet^ 


-11.55 ±0.04 


-11.72 ±0.08 


Si IV doublet* 


-12.20 ±0.71 


-12.39±0.05 


N V doublet* 


-11.79 ±0.04 


-12.21 ±0.22 


LycK* 


-11.32 ±0.04 


-11.23 ±0.04 



^Parameters of a simple power law fit, of the 
form oo X (A/1750 A)~", where ao is the contin- 
Tium flux density at 1750A (in units of 10~^^ erg 
s^^ cm~^ and a is the power law index. 

*The logarithm of the total line flux 
(ergss~^ cm~^) is shown for the indicated 
doublet or line. 



Table 8. Swift XRT Observations 



Obs ID 


UT Date 


Exp Time 


Count Rate* 






(s) 


(ctss-^) 


00036620001 


2007-05-12 


6670 


0.56 ± 0.09 


00036620002 


2007-07-05 


2355 


0.70 ± 0.02 


00036620003 


2007-09-19 


2213 


0.22 ±0.01 


00036620004 


2007-09-20 


2824 


0.25 ±0.01 


00037592001 


2009-06-14 


5149 


0.77 ±0.01 



*Count rate is given in the 0.3-10 keV band. 
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Table 9. Best-fit Parameters for the X-ray Observations 



Obs. Date Akr'' Pt F(0.3-10keV)t x^/dof 



2007-05-12 (XRT) 


1 07+1.46 
-L-'J' -1.34 


1 07+0.07 


1.98 ±0.07 


169.7/144 


2007-07-05 (XRT) 


1.79* 


1 9+0.05 
^•-"--^-0.06 


2.50 ±0.12 


63.7/73 


2007-09-19 (XRT) 


9 00+0.99 
^■^^-1.14 


-'■•'3 "3-0.16 


0.88 ±0.10 


14.6/21 


2007-09-20 (XRT) 


^■'-'"-0.99 


^^■"^"-0.12 


-■^■^^--0.09 


41.9/32 


2009-06-14 (XRT) 


1.95* 


2 in+o-03 

^•-'-'-'-0.04 


9 71+0.08 
^- ' ^-0.10 


131.3/150 


2009-12-13 (pn) 


14 q+i-34 


2.08 ± 0.03 


3.34 ± 0.06 


391.9/339 



^Parameters are shown from a blackbody + power-law fit. Tlic blackbody 
temperature was fixed at 0.10 keV and the flux normalization is recorded as 
AkT (in units of 10~^ photons keV"-*^ cm~^ s~^ at 1 keV). Parameters of the 
power-law (using the pegpwrlw model) include the photon index, T where 
F{E) oc E~^^~^\ and F(0.3-10keV), the flux in the power-law component 
in the 0.3-10 keV band in units of 10~^^ ergs~^ cm~^. 

■'■Upper limit on the indicated parameter. 



Table 10. Measured Flux from the UVOT/OM Observations 



Obs. Date 


ut 


uvwit 


UVM2t 


UVW2t 


2007-05-12 (UVOT) 


1.75 ± 0.07 








2007-07-05 (UVOT) 


1.64 ±0.07 


2.38 ±0.11 






2009-06-14 (UVOT) 


2.45 ±0.10 


3.64 ±0.17 


5.00 ±0.15 


6.17 ±0.27 


2009-12-13 (OM) 




2.88 


3.74 





^Fluxes are in units of 10 ergs ^ cm" "2 A-^ For the UVOT, the central 
wavelength of the filters corresponds to U at 350lA, UVWl at 2634A, UVM2 
at 2231A, and UVW2 at 2030A. For the OM, the filters correspond to UVWl 
at 29IOA and UVM2 at 2310A. 
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Fig. 1. — Spectra of the UV observations from lUE, GHRS, and COS are shown, labehng the 
locations of emission lines frequently seen in AGN. The spectra are corrected for Galactic extinction. 
Additional absorption features in the spectrum are mostly associated with absorption within our 
galaxy. The most prominent intrinsic AGN features are labelled. 
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Fig. 2. — Emission line fits to the COS ERO spectra are shown. The C IV doublet is shown at top, 
the Si IV doublet, blended with O IV], in the middle, and Lya + N V in the bottom panel. The 
flux (black) and errors (green), are normalized by the best-fit power law continuum. The dashed 
lines show individual Gaussian components, while a solid red line shows the combined best-fit 
emission line model. Galactic ISM lines are labeled; the most prominent absorption features were 
masked for fitting. Additionally, the COS spectra were smoothed for illustrative purposes using 
moving- window averaging of four pixels. 
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Fig. 3.— Same as Figure El but for the COS GTO data. 
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Fig. 4. — Plotted are the Lya emission line regions for the 2009 COS ERO (green) and GTO (blue) 
observations along with the 1997 GHRS spectrum (black). The ERO observation had the highest 
flux level. The inset shows a zoom in around the potential outflowing Lya absorption feature that 
was detected in the GHRS spectrum, as well as the four 2000-2001 FUSE observations. Clearly, 
the feature is weaker in the COS observations (by a factor of ~ 5). Additional absorption features 
in the spectrum include Galactic I SM features ( r ed) an d HVCs /additional w eak intrinsic absorbers 
(magenta), which are detailed in iPenton et al.l (|2000l ). I Collins et al.l (|2003l ). and this paper. See 
Figure El and the text for more details on these features. 
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Fig. 5. — Comparison of the normalized GHRS (red) and COS ERO+GTO data (black) is shown. 
Lya absorption features noted by Penton et al. (2000) are marked with filled stars and open stars 
mark three additional features noted here (two absorbers, Aobs = 1223.5 A and 1224.2 A, are not 
shown). Absorption arising from interstellar N V and S II lines are marked. Several features appear 
unchanged between observation epochs while several others exhibit marked changes. See text for 
discussion of individual systems. 
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Fig. 6. — Plotted is a comparison of the spectra and best-fit model for the COS ERO (green) 
and lUE (magenta) spectra in C IV and the Lya + N V regions. The best-fit emission line fits 
are shown in black. The COS spectra are smoothed by moving-window averaging of four pixels. 
Low-resolution in the lUE spectrum coupled with strong geo-coronal Lya emission make it difficult 
to compare the line flux in Lya -|- N V doublet with the higher resolution COS spectra. 
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Fig. 7. — The UV light curve for Mrk 817 is plotted, showing changes in luminosity (in ergss^^) for 
the emission lines (Lya + N V doublet, C IV doublet, and the Si IV doublet, blended with O IV]) in 
the top panel and the continuum flux density (the normalization value oq from the power law fits) 
in the bottom panel. The horizontal lines show the averages for the lUE and COS observations 
in each of the values. The GHRS measurements are estimates based on spectral fitting in the 
Lya range, scaling to the broader bandpass fits of the COS ERO observations as explained in the 
text. We find the lUE observations to show lower continuum and line luminosities than the COS 
and GHRS observations. However, there are no drastic changes over the observations, which span 
almost 30 years. 
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Fig. 8. — The XMM-Newton pn spectrum of Mrk 817 is shown in the observed 0.4-2 keV band, 
re-binned to a signal-to-noise ratio of 10. This X-ray observation was taken within two weeks of the 
COS GTO observation. The bottom panel shows the ratio of the data/model, where the model is 
a blackbody + power-law model. The low-energy spectrum exhibits no signs of an outflow, which 
is what we expected based on the lack of a strong detection in the UV COS observations. The 
location of the O VII and O VIII absorption edges are indicated on the plot. 
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Fig. 9. — Plotted is the ratio of the data/model for the blackbody + power-law model fit to 
the December 2009 XMM-Newton pn spectrum (black, re-binned to a signal-to-noise of 20 for 
illustrative purposes) and the five Swift XRT spectra (re-binned to a signal-to-noise from 5 to 10, 
depending on the exposure time of the observation). The model applied to the spectra is the best-fit 
model to the pn spectrum. It is clear that the source varied both in flux and spectral shape between 
the X-ray observations. The lowest flux points (both in blue) correspond to the Swift XRT spectra 
taken a day apart in August 2007. The additional Swift spectra correspond to the May 2007 (red), 
July 2007 (green), and June 2009 (magenta) observations. 
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Fig. 10. — The photon index versus 0.3-10 keV luminosity for the X-ray observations of Mrk 817 
is plotted, which span a range of approximately 2.5 years. There is a clear correlation between 
the luminosity of the source and the photon index, with steeper photon indices corresponding to 
higher luminosities. We plot the result of a linear regression fit, which shows a strong {B? = 0.94) 
correlation between F and Lx- 
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Fig. 11. — The X-ray light curve of Mrk 817 is plotted, including the 0.1-2.4 kcV luminosity mea- 
sured from two ROSAT PSPC observations, five Swift XRT observations, and one XMM-Newton 
observation. The X-ray luminosity of Mrk 817 is clearly variable, with the highest luminosity 
observation 39 times that of the lowest luminosity observation. 
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Fig. 12.— Plotted are the UV-optical fluxes obtained from the Swift UVOT and XMM-Newton OM 
observations that were simultaneous to the X-ray observations. Additionally, the best-fit emission 
line (solid lines, shown in the rest wavelength of Mrk 817) and continuum fits from our analyses of 
the COS spectra are shown (dashed lines), extended into the optical band. The UVOT and OM 
fluxes are at a similar level to the values predicted by the COS observations, particularly for the 
OM observations which were taken within two weeks of the COS GTO spectra and are within 10% 
of the values predicted by the COS GTO continuum. 
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Fig. 13. — Plotted are the spectral energy distributions from XMM-Newton and Swift observations 
with associated optical/UV data from OM/UVOT. A model consisting of an accretion disk model 
with an absorbed power-law component was fit to each of the individual observations. The opti- 
cal/UV portion of the SED does not change significantly, while there is a noticeable difference in 
the X-ray portion. 



